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Abstract. How the solar wind is accelerated to its supersonic speed is intimately 
related to how it is heated. Mechanisms based on ion-cyclotron resonance have been 
successful in explaining a large number of observations, those concerning the signifi- 
cant ion temperature anisotropy above coronal holes in particular. However, they suffer 
from the inconsistency with turbulence theory which says that the turbulent cascade 
in a low-beta medium like the solar corona should proceed in the perpendicular rather 
than the parallel direction, meaning that there is little energy in the ion gyro-frequency 
range for ions to absorb via ion-cyclotron resonance. Recently a mechanism based on 
the interaction between the solar wind particles and the anisotropic turbulence has been 
proposed, where the perpendicular proton energy addition is via the stochastic heating 
(Chandran et al. 201 1). We extend this promising mechanism by properly accounting 
for the effect of proton temperature anisotropy on the propagation of Alfven waves, 
for the radiative losses of electron energy, and for the field line curvature that naturally 
accompanies solar winds in the corona. While this mechanism was shown in previous 
studies to apply to the polar fast solar wind, we demonstrate here for the first time that 
it applies also to the slow wind flowing along field lines bordering streamer helmets. 



1. Modeling the solar wind: why temperature anisotropy? 



How the solar wind is accelerated is closely related to how it is heated (e.g.. lCranmer 
20091) . For the nascent fast wind above coronal holes, much evidence exists that ions 



are hotter than electrons, and their perpendi cular temperatu res a re higher than the par 



allel one (T 1 - > r") (see data compiled in ICranmer 2002 . also Cranmer et al. I l2008h 



Interestingly, in the near-Sun region the slow wind flanking streamer h elmets exhibits 



a sim i lar tendency as ev idenced by measurements with SOHO/UVCS dStrachan et al. 



20021 : iFrazin et al.1 120031) . Given that the solar winds receive most acceleration close 



to the Sun, these observations point to the possibility that the fast and slow winds may 
share a common heating, and hence acceleration, mechanism. Some other factors, most 
notably the distinct flow tube geometry, then determine the h eating profile and hence 
lead to significantly different terminal speeds of the two winds ( Wang & Sheeleylll990l : 
ICranmer et"al"]|2007h . 

The temperature measurements, especially the inferred significant ion tempera- 
ture anisotropy, lead naturally to the suggestion that the nascent solar winds are heated 
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by ion-cyclotron resonance (e.g.. iHollweg & Isenbergl [20021) . The needed high fre- 
quency ion-cyclotron waves (ICWs) may be generated either by a turbulent parallel 
cascade from low-frequency waves emitted by the Sun, or directly by small-scale mag- 
netic reconnection events at chromospheric network. While the latter remains safe, 
the former has bee n challenged rec ently as both MHD tubulence theory and simu- 
lation studies (e.g.. ICho et all 120031) suggest that in a low-/? solar corona, the turbu- 
lence cascade prefers the J_-direction, hence generating modes with high k ± instead 
of ICWs with high k\\ and high frequency. However, these high-& ± modes, presum- 
ably kinetic Alfven waves (KAWs), are most readily dissipated by electrons in a low-/? 
plasma, if the only dissipation channel is via the Landau resonance. The consequent 
preferential electron heating would violate available measurements. A remedy comes 
only recently, where KAWs with sufficient amplitudes were demonstrated to primar- 
ily heat protons, preferen tially in the ^-direction, by rendering the ion orbits stochas- 
tic dChandran et al.ll2010h. This nonline ar mechanism, termed "stochastic heating", was 



further paramete rized (Chandran 2010h and eventually incorporated into a two-fluid so- 
lar wind model dChandran et al.l 201lh (hereafter Cll) where the anistropic turbulence 
was shown to be capable of producing polar fast solar wind solutions that agree with a 
substantial number of observations. 

Given that the fast and slow winds may share a common heating mechanism, and 
that CI 1 is restricted to the polar fast wind, we wish to extend CI 1 by asking whether 
the mechanism therein also applies to the equatorial slow wind as well. To address 
this, we construct slow wind models that are based on the anisotropic turbulence, that 
disinguish between electrons and protons, and that account for the proton temperature 
anisotropy. Inherited from Cll are: 1) The turbulence is driven by reflections of out- 
ward (z + ) waves off the gradients of the solar wind parameters, 2) Outward waves z + 
dominate inward ones z~, 3) z~ waves are cascaded sufficiently fast compared with their 
wave periods, 4) the cascade proceeds in the inertial range along a "critical balance" 
path until it reaches a J_-scale comparable to the ion gyroradius, where the generated 
KAWs are damped via both Landau damping (yielding electron heating, and proton 
H-heating), and stochastic heating (yielding proton -L-heating), 5) if the wave energy is 
cascaded to smaller still X-scales, only electrons receive this energy. However, com- 
pared with Cll, we: 1) properly incorporate in the gross dissipation rate the field line 
curvature that accompanies field lines bordering streamer helmets, 2) incorporate the 
radiative loss in the electron energy equation, which is important in the energy balance 
in the upper transition region, 3) improve the treatment of wave propagation by includ- 
ing the effects of temperature anisotropy, which was missing in the original treatment 
in Cll. For simplicity, instead of a more self-consistent, collisionless treatment of the 
proton heat fluxes based on the Landau fluid approach as was adopted by Cll, we as- 
sume that the fluxes can be described by a modified Spitzer form. In what follows, we 
describe the improved model in section |2l present a constructed slow wind model in 
section |3j and summarize our effort in section |4] 



2. Model desciption and method of solution 

The fluid pa rt of our model starts with the standard 16-moment transport equations (see 
appendix in iLi & l1|2009|) . When axial symmetry is assumed and solar rotation ne- 
glected, in a steady state the vector equations may be decomposed into a force balance 
condition across the meridional magnetic field lines and a set of transport equations 



Anisotropic-turbulence-heated solar wind model 



3 



along them. For simplicity we r eplace the former by prescri bing a meridional magnetic 



field configuration presented in lBanaszkiewicz et al.l (1 19981) . which is representative of 



a minimum corona. The transport equations along curved field lines then read 

inva)' = 0, (1) 
k B r / || x V k B (Tl-T£)a' /GMnV F 

vv' + 



[n{4 + Te) J + _p _J_ =0 , (2) 

y L v /J arrip \ r I nm p 

v(T e ) + — [aK M T e ' T e ) 

a nk B a v ' 

+ 2v pe (T e - T p ) - ^-{Q e - L md ) = 0, (3) 
nk B 

ii I, 3(7 - 1)0" 

+ 2v pe (Tl - T e ) + 2v pp {Tl - T£) - P = 0, (4) 

. T±v 3(y-l)vrta / y - 1 r , y v 

3(r- 1)<2£ 

+ iv^cr^ - r e ) + Vpp (T£ - rh — — '- = o. (5) 



Here n is the number density and v the speed, a <x 1 /B is the tube cross-sectional area 
with B being the magnetic field strength. An arbitrary point along a field line is charac- 
terized by both (r, G) (its heliocentric distance and colatitude) and arclength /. The prime 
' denotes the differentiation with respect to /. Moreover, k B and m p are the Boltzmann 
constant and proton mass. The gravitational constant is denoted by G, and solar mass 
by M Q . The electron, proton parallel and perpendicular temperatures are denoted by T e , 
T p and T p , respectively. The mean proton temperature is given by T p = (T p + 2T p )/3. 
Furthermore, y = 5/3 is the adiabatic index. The Coulomb collision rates v pe and v pp 
are evaluated by using a Coulomb logarithm of 23. For simplicity, the Spitzer law is 
assumed for the electron heat flux, K e o = 7.8 x 1(T 7 (cgs units will be used throughout). 
A similar form is assumed for the proton ones, but the coefficient is arbitrarily reduced 
relative to the Spitzer value tc pc = 3.2 x 10~ 8 by a spatially varying factor, namely, 
kpo/Kp,c is 1 for r < 3 R n , rampe d linearly with r to 0.02 at r = 10 R Q and kept so from 
there on (see e.g., iLi et al.|[l999l) . B esides, L raf i represen ts the radiative losses and we 
adopt the standard parametrization by lRosner et all d 19781) for an optically thin medium. 
2e> Qp an d Q p are the heating rates, while F represents the volumetric force density 
the solar wind receives. Part of the electron heating is assumed to come from some 
ad hoc process which operates close to the base, <2 ei basai = QeO exp(-/// f /). This basal 
heating is found to be important to improve the agreement of the model results with 
observations. On the other hand, the majority of the electron heating <2e, W av, as well as 
the proton heating rates, are due entirely to the turbulent dissipation of low-frequency 
Alfven waves. 

The wave part of the model is a WKB evolution equation with dissipation terms, 

+ V F = -£ wav . (6) 



4 



Li & Habbal 





10 5 




1Q -12 




10 -13 


w 






io- 14 








10 -15 








10 -16 


io- 17 




1.000 


O 


0.100 




0.010 



0.001 



1 



eI 


r Q P 7<?tot 


(d) 1 




\ Qe/^tot 


0^0^ 



10 



100 



-(R s ) 



Figure 1 . Radial distribution of solar wind parameters derived from an anistropic- 
turbulence-heated 2-fluid model with proton temperature anisotropy. (a): number 
density n and speed v. (b): electron, proton parallel and perpendicular temperatures 
T e , Tp, and T£. (c): gross turbulence dissipation rate per particle Q^/n and the 
geometrical factor ln(r sin On ' 4 ). (d): fractions of gwav that goes to heating electrons 
and protons. In (a) and (b), the error bars and open boxes represent some remote- 
sensing measurements representative of solar minimum conditions (please see text). 
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Here only the outward waves contribute to the wave energy flux density F w and wave 



force F, whose expressions are given by Eqs.(18) and (16) in lLi et aLl(ll999h . where the 
correction due to finite temperature anisotropy is included. On the other hand, the gross 
dissipation rate does take into account the inward waves, and the interaction between 
outward and inward ones yields 

Gwav = 4c dPw (v + v A ) | [in (r sin #« 1/4 )]'| , (7) 

where p w is the wave pressure, va the Alfven speed, and is a dimensionless constant. 
Note that the field line shape appears here by contributing further reflection via the 
gradient of the geometrical factor sin#. To apportion Q wav among <2e,wav and (Jp ± , 
Eqs.(44) to (56) in Cll are used. It suffices to mention here that a correlation scale 
L i oc 1 / V# appears not directly in Q wav , as in ICW-based models ( Hollweg & Isenbergl 



2002), but in the way <2 wav is apportioned. 

To proceed, the field line flanking some streamer helmet is characterized by the 
B-distribution and its shape. T hese are imposed by requiri ng that the field line reach 



colatitude at 1 AU in the iBanaszkiewicz et all (1 19981) model with their fiducial 



parameters (the DQCS model as presented in Fig. 3 therein), and that B at 1 AU be 
3.5 y. The wind is determined by the following set of parameters: the base density 
n e , temperature T Q , and wave amplitude 6v Q , as well as L ± at the base L©. We choose 
n Q = 2.5 x 10 8 cm" 3 , T Q = 1 x 10 5 K, Sv Q = 35 km s~\ and L Q = 100 km. Once 
these are prescribed, we first cast equations CD to © in a time-dependent form and then 
evolve them from an arbitrary initial state until a steady state is found, which is then 
examined in some detail in next section. 



3. Constructed slow solar wind model 

Figure Q] presents the radial distribution between 1 R Q and 1 AU of a number of solar 
wind parameters. In Fig{T^, the compute d number dens i ty n is compared with the elec- 
tron density measurements presented in IStrachan et al.l (l2002h (their fig.3c), while the 
speed profile v is compared with the range of wind speeds derived by tracking a c ollec- 



tion o f small inhomogeneities (the blobs) in images obtained with SOHO/LASCO (IWang et al 
l2000h . In addition to obtaining reasonable values of speed v = 342 km s" 1 and flux den- 
sity nv = 3.5 x 10 s cm -2 s" 1 at 1 AU, the model results agree reasonaly well with the 
observations close to the Sun. The temperatue profiles, presented in FigflJ), are found 
to be difficult to reproduce th e measured H I ones obtained with SOH O/UVCS and 



repre sented by the open boxes ( Strachan et al.ll2002l) as well as error bars dFrazin et al 



20031) . Despite this, the resultant pressure gradient force produces a decent slow wind 



solution. It is interesting to see that, the gross turbulence dissipation rate <2 wav exhibits 
multiple zeroes (FigQJ, black curve), as a result of the geometrical factor ln(r sin dn 1 ^ 4 ) 
(FigHfc, red curve) attaining its local extrema. This is understandable given that <2 wav is 
proportional to the gradient of this factor (Eq|7]). Also noteworthy is that, among the dis- 
sipation channels of the wave energy, the Landau damping of KAWs does not make any 
significant contribution (FigOJl). When the electron heating dominates (r < 1.6 R Q ), 
the dissipation is due to the processes that take place at perpendicular scales shorter 
than the proton gyroradius. On the other hand, when the proton heating dominates 
(r > 1.6 Rq), the dissipation is almost entirely due to stochastic heating of KAWs, 
evidenced by the fact that protons literally receive no parallel heating. At distances 
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r > 122 Rq where Tp-/T p is substantially smaller than unity, the solution becomes 
firehose unstable. However, this should have little effect on the gross wind parameters 
given that the wind itself has already been fully accelerated. 



4. Summary 



So far the rather promising mechanism, originated by iMatthaeus et aL ( 19991) a nd fur - 
ther developed by, t o name but a few. IVerdini et all (120051): IVerdini & Velli! (I2007Q: 

Verdini et al.l(l20100 : ICranmer et all (120071) : ICranmer & van Ballegooijenl ( 2*0121) : IChandran et al 
(1201 II) . is almost entirely devoted to examining the polar fast solar wind. Does it also 
apply to the slow wind that flows along curved magnetic field lines bordering streamer 
helmets? Our preliminary results suggest that this is indeed the case. 
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